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Abstract. We carried out AC magnetic susceptibility measurements and muon 
spin relaxation spectroscopy on the cubic double perovskite Ba 2 YMo06, down to 50 



I 

mK. Below ~ 1 K the muon relaxation is typical of a magnetic insulator with a 
spin-liquid type ground state, i.e. without broken symmetries or frozen moments. 
i However, the AC susceptibility revealed a dilute-spin-glass like transition below 

~ 1 K. Antiferromagnetically coupled Mo 5+ 4c? 1 electrons in triply degenerate t 2g 
orbitals are in this material arranged in a geometrically frustrated fee lattice. Bulk 
magnetic susceptibility data has previously been interpreted in terms of a freezing to a 
00 heterogeneous state with non-magnetic sites where 4c? 1 electrons have paired in spin- 

singlets dimers, and residual unpaired Mo 5+ 4c? 1 electrons. Based on the magnetic 
heat capacity data it has been suggested that this heterogeneity is the result of kinetic 
constraints intrinsic to the physics of the pure system (possibly due to topological 
overprotection), leading to a self-induced glass of valence bonds between neighbouring 
4c? 1 electrons. The ^.SR relaxation unambiguously points to a static heterogeneous 
state with a static arrangement of unpaired electrons isolated by spin-singlet (valence 
bond) dimers between the majority of Mo 5+ 4c? electrons. The AC susceptibility data 
indicate that the residual magnetic moments freeze into a dilute-spin-glass-like state. 
This is in apparent contradiction with the muon-spin decoupling at 50 mK in fields up 
to 200 mT, which indicates that, remarkably, the time scale of the field fluctuations 
from the residual moments is ~ 5 ns. Comparable behaviour has been observed in other 
geometrically frustrated magnets with spin-liquid-like behaviour and the implications 
of our observations on Ba 2 YMo06 are discussed in this context. 
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1. Introduction 

Geometrical frustration of the exchange interactions in antiferromagnetic (Mott) 
insulators with magnetic topologies based on triangles and tetrahedra can fully 
suppress the spontaneous symmetry breaking to the common long-range ordered 
antiferromagnetic (Neel) state [U E]. Quantum-ordered spin-liquid states are widely 
expected to occur in such system. These could be gapless, with emergent Fermionic 
quasiparticles called spinons with S = 1/2 and no charge or gapped (topologically- 
ordered) [31 H] . A number of materials with a spin-liquid-like absence of frozen moments 
have now been identified [SI El El El E] • These are without exception gapless, often still 
display some weak hysteretic behaviour below 2 K, and have only short-ranged dynamic 
magnetic correlations [TU1 [HI IB] which is contrary to theoretical expectations of gapless 
spin liquids. 

While these short-ranged correlations could in some cases be due to disorder [101 18] , 
another plausible explanation is the presence of kinetic constraints towards freezing 
into the (quantum) ground state [121 [131 El EE]- It has been shown that even where 
quantum fluctuations are so strong that no local order parameters exist, as in some 
geometrically frustrated quantum magnets, these quantum fluctuations are not always 
effective for equilibration to its true quantum ground state [121 H31 El US]- When 
the width of the tunnelling barriers becomes comparable to the system size, out-of- 
equilibrium quantum-glassy phases can stabilize [121 EES E] • Quantum spin systems with 
strictly local interactions, with local degrees of freedom coupled to a heat bath (Mott 
insulators), also studied in the context of fault-tolerant quantum computation [To] , are 
ideal model systems to study this quantum glassiness [12]. A particularly intriguing 
and far-reaching possibility are that there are kinetic constraints due to "topological 
overprotection" [T21 EH] . The topological order [17] in a gapped quantum (spin liquid) 
state is known to lend the ground state a certain rigidity and "protect" the quantum 
wavefunction against decoherence and material imperfections [H [16]. An infinite 
number of local operations would be required to mess up the macroscopic entanglement 
of this quantum ordering. In Mott insulators this concept can be turned on its head; 
here the physics is local [19J and hence the only route to access any topologically- 
ordered states in the first place, is via an infinite number of local operations. This has 
been termed "topological overprotection" [T21 EH] . Within this context we studied the 
cubic fee lattice antiferromagnet Ba2YMo06, which has previously been suggested [20] 
to freeze into a valence bond glass (VBG) state [21]. This is in line with theoretical 
predictions for a quantum magnet with an fee lattice predicting a "strong glass" state due 
to topological overprotection [121 HI]- However, the Hamiltonian of this cubic perovskite 
with unquenched orbital degrees of freedom is significantly more complex than that of 
the Heisenberg Hamiltonian, as was recently pointed out |22j . 

Ba2YMoO"6 is a double perovskite (of general formula A2BB'Oe) where the two 
B-site cations have assumed rocksalt order [23] (Fig. la). Magic-angle spinning NMR 
was reported to point to 3 % cation disorder [24J but in neutron diffraction data of our 




Figure 1. (Color online) (a) The Ba2YMo06 double perovskite structure, showing 
the rocksalt ordering of MoOg octahedra (blue) and YOg octahedra (yellow) with Ba 2+ 
(green) at the perovskite A site, (b) The network of edge-sharing tetrahedra of the 
fee antiferromagnet. (c) An illustration of the valence bond glass state, with the non- 
magnetic spin singlet dimers across random edges (yellow) and a single unpaired spin 
(red), (d) The muon stopping site (red) on an edge of the YO@ octahedron, surrounded 
by 6 Mo0 6 octahedra. 

sample the B-site ordering of the Mo/Y cations refines to 100(1)% 15, as expected 
for the large difference in ionic radii of the Y 3+ and Mo 5+ of 0.29 A. The Mo 5+ 
cations have a single unpaired electron in an octahedral crystal field and are separated 
from each other by 90° Mo-O-O-Mo bonds (J\) along the edges of corner-sharing 
tetrahedra, and 180° Mo-O-Y-O-Mo bonds (J 2 ) (Fig. lb). The 12 near-neighbour (Ji) 
bonds and 6 further neighbour bonds per site lead to a strong net antiferromagnetic 
interaction, with a Weiss temperature of -143(5) K. Ba2YMo06 remains cubic (space 
group FM3m) down to at least 2 K which means the single 4d t 2g electron is also 
orbitally degenerate. Due to relatively strong spin-orbit coupling in the 4d shell this 
degeneracy is expected to be partially lifted, leading to J = 3/2 quadruplets [22J. 
It is noteworthy that of the related compounds where the Y 3+ is replaced for a 
trivalent lanthanide ion only in the case of Nd 3+ and Sm 3+ a Jahn- Teller distortion 
is observed to a tetragonal unit cell, and in both cases this distortion is accompanied 
by antiferromagnetic ordering |23j. The actual Hamiltonian that applies to Ba2YMoC>6 
is not reproduced here for brevity. It contains antiferromagnetic and ferromagnetic 
exchange interactions that depend on the orientations of the occupied orbital states 
within the degenerate ti g manifold, a quadrupolar interaction between neighbouring 
orbitals and a spin-orbit coupling term which can be projected-out by constraining the 
system to a low-energy subspace with J = 3/2 moments [22J. A (naive) estimate of the 
orbital-dependent 15 nearest-neighbour antiferromagnetic exchange interaction strength 
Ji, assuming further-neighbour interactions can be neglected, yields | J\\ ~ 72 K. 

Remarkably, the magnetic/electronic heat capacity of the single 4c? t2 9 electron 
at the Mo site evidenced a pseudogap below ~ 50 K, gradually locking up most of 
the magnetic entropy of the J = 3/2 moments [20J. This "strong glass" -like freezing 
does however not translate to a significantly increased muon relaxation down to 2 K, 
indicating that most of the Mo 5+ 4d electrons pair up with near-neighbour spins to 
form non- magnetic spin-singlet- or valence bond dimers [20J (Fig. lc). A second 
Curie-Weiss regime is observed in the DC susceptibility below 25 K, with a Weiss 
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temperature of —2.6(2) K. The paramagnetic susceptibility between 1.9 and 7 K was 
found to be equivalent to that from 8(1)% of the Mo 5+ magnetic moments as measured 
at high temperature (see supplementary information). These residual, apparently 
weakly-coupled spins (sometimes referred to as "orphan spins") are often observed in 
geometrically-frustrated compounds [2S]- In some cases this has been suggested to 
be spins left-out from spin-singlet valence-bond dimerizations between the majority 
of spins, leaving a small fraction of magnetically-isolated and thus weakly-interacting 
unpaired-spin moment |26| [TO] (Fig. lc). Hence, "dangling spins" might in this case be 
a better term. The aim of the present study was to characterise the low-temperature 
spin-orbital dynamics to confirm the VBG proposal. 

2. Results 

The polycrystalline sample was synthesised using a solid state synthesis method as 
described in [20]. AC magnetic susceptibility measurements were carried out using 
an inductive susceptometer from CMR-direct in an Oxford Instruments dilution fridge 
(down to 50 mK) and in a 3 He fridge (down to 300 mK). The measurements used an 
AC field of 0.02 mT in a frequency range from 54.5 to 8900 Hz. The intensity of the AC 
susceptibility signal at all frequencies was assumed to match the low-temperature DC 
susceptibility obtained using a SQUID magnetometer [see supplementary information]. 
As shown in Fig. 2a,b frequency-dependent cusps are observed around 600 mK in 
both the dispersive (%') and dissipative (%") parts of the dynamic susceptibility. The 
frequency dependence of T max (Fig. 2c) exhibits Arrhenius-activated behaviour (the 
straight line in Fig. 2c), with an activation energy E a = 39(2) K and a trial frequency 
f ~ 1 x 10 29 Hz. Such a frequency dependence in the AC susceptibility is a defining 
feature of a strong glass as well as spin glasses [27] . Figure 2d shows the effect of an 
additional DC field on the AC susceptibility at the spin-glass transition, revealing only 
a modest suppression of the AC susceptibility in fields below 250 mT. The AC magnetic 
susceptibility is a macroscopic property. We complemented this data with muon spin 
relaxation spectroscopy which is a local probe of the spin dynamics in the material, 
using the MUSR instrument at the ISIS facility, UK. Figure 3a shows the average time 
dependence of the spin polarisation of muons implanted in a 3 g powder sample, for a 
selection of temperatures between 50 mK and 120 K and in zero external field. Figure 3b 
shows the field dependence of the muon relaxation at 50 mK, for (longitudinal) magnetic 
fields applied along the initial muon spin polarisation direction. 

At the time the muons are implanted in the sample (t — 0), all muon spins point 
along the forward (z) direction. The muon spins will then precess due to the local fields 
at the muon sites. In a material with uniform magnetic order all muons will precess 
with the same frequency, which translates to an oscillation of the muon polarisation 
P z {t)- In the case of static disordered moments a single damped oscillation is observed 
equilibrating to l/3rd of the initial polarisation, corresponding to on average l/3rd 
of the randomly orientated magnetic moments that are aligned with the initial muon 
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Figure 2. (Color online) AC magnetic susceptibility measured in a dilution 
refrigerator. (a) The frequency dependence of dispersive (real) part of the AC 
susceptibility at low temperatures matched to the DC susceptibility. (b) The 
dissipative (imaginary) part of the AC susceptibility near the spin-glass transition, 
(c) The Arrhenius- activated frequency dependence of the temperature with maximum 
AC susceptibility (T max ) (see text) (d) The effect of DC fields up to 1 T on the AC 
susceptibility (at 1029 Hz) near the transition temperature. 

spin polarisation and hence do not cause muon precession or depolarisation. This type 
of relaxation is observed, for example, below the spin-glass transition in the related 
compound S^MgReOe with Re 6+ S = 1/2 [28]. The muon relaxation in the presence of 
static but random fields can be suppressed by the application of a longitudinal magnetic 
field stronger than the local fields at the muon site. This enabled the determination of 
the contribution to the muon relaxation from nuclear spins P n uc{t) & t the Y and at a 
fraction of the Ba and Mo sites in Ba2YMoC>6, which below ~ 10 K are static at the time 
scale of the muon experiment. At 1.3 K, just above the spin-glass transition observed in 
the AC susceptibility data, a 5 mT field was sufficient to decouple the muon relaxation 
for the small nuclear fields (typically up to 1 mT ), yielding the nuclear relaxation factor 
Pnuc{t) as shown in Fig. 3b. The observed nuclear relaxation is as expected for the 
muon stopping site on the edges of the Y06 octahedra [29] (Fig. 3d, see supplementary 
information for further details). No further change in the muon relaxation was observed 
with the application of longitudinal fields up to 200 mT, confirming that the remaining 
muon relaxation was due to larger (non-frozen) electronic moments. Well into the 
paramagnetic state the electron-spin fluctuations are usually too fast to be picked-up 
by muons leading to negligible muon relaxation at 120 K (Fig. 3a). As indicated 
by the heat capacity data, the entropy from most of the J = 3/2 moments is locked 
up below ~ 40 K. The very slow relaxation at 2 K thus confirms that the majority 
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Figure 3. (Color online) Muon spin relaxation data obtained at MUSR, ISIS, (a) The 
temperature dependence of the zero- field muon relaxation (P z (t)) between 50 mK and 
120 K. (b) The field dependence of the muon relaxation at base temperature (50 mK) 
and the nuclear relaxation factor which was confirmed to be constant below 1.3 (See 
text and supplementary material), (c) The fit parameters for fits of Eq. 1 to the zero- 
field data as a function of temperature, taking into account the weak nuclear relaxation, 
(d) The fit parameters for fits of Eq. 1 to the 50 mK muon relaxation in fields up to 
200 mT. The fluctuation rate was determined from the field dependence of Ai using 
Redfields equation (solid green line) (see text and supplementary information). 



of these magnetic moments must have paired into non-magnetic spin-singlet or valece 
bond dimers [20] ■ Across the spin-glass transition the muon relaxation rate increases 
sharply but even at 400 mK, below which no further increase in the muon relaxation 
(depolarization) is observed, the relaxation rate remains remarkably low, even for a 
system with dilute electronic spins. Amongst antiferromagnetic insulators without a 
spin-gap, only in the spin-liquid compound herbertsmithite is a slower muon relaxation 
observed at 50 mK [6J. If this weak relaxation were due to static random moments, 
then the fields from these moments at the muon site would have to be as small as 
~ 0.6 mT, so that the muon relaxation should be fully decoupled in a 5 mT field. As is 
clear from Fig. 3b, a 5 mT field only leads to a small reduction of the relaxation, the 
difference with the zero-field relaxation being given approximately by P nuc {t). A clear 
muon relaxation signal persists in fields up to 200 mT, consistent with the presence of 
fluctuating electronic magnetic moments. 

The muon relaxation P z {t) is determined by the magnitude, variation and 
fluctuation of local fields inside a sample, and carefully analysing the shape of the 
muon relaxation curve can allow to distinguish different scenarios. Here we found that 
the best fitting scenario is with two different muon environments, designated El and 
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E2, with an exponential relaxation in El (arising from spin fluctuations with a rate Ai 
that will be estimated later), and a phenomenological modified-exponential decay with 



where P nuc (t) = 1 for longitudinal fields ^ 5 mT, / is the fraction of muons in El and 
the constant background b = in zero-field. The results of the fits to the data down 
to 50 niK and in fields up to 200 mT are shown in Fig. 3c and d, respectively. The 
background contribution b in E2 is constrained by the condition that with the fields 
up to 200 mT applied here the fraction of muons in El and E2 should be constant, as 
implied by the AC susceptibility data in Fig. 2d. Just above the spin glass transition, at 
1.3 K, the stretching exponent /3 = 1 and Ai = 0. This implies that a major fraction / 
of the muons is located in a non-magnetic environment (El), while the remaining muons 
(~ 30% at 1.3 K) undergo an exponential relaxation due to surrounding fluctuating spins 
(E2). On cooling through the spin-glass transition the depolarisation of the muons in 
E2 changes sharply from exponential to Gaussian, as previously observed in the S = 3/2 
kagome antiferromagnet SrCrgGa^ig [26], and the fraction of muons in E2 increases 
to 40(5)%. At the lowest temperatures accessed the muons in El also exhibit a weak 
exponential relaxation (Fig. 3c). At base temperature El and E2 correspond to muons 
surrounded by respectively and 1 unpaired electron on the 6 Mo 5+ ions surrounding 
the muon site. Non-magnetic (El) and magnetic (E2) environments were also evident 
from 89 Y NMR [24j with a ratio of approximately 50:50, comparable to our results (Fig. 
3c). This confirms that 1) the muon stopping site is near the Y 3+ and 2) that the muon 
does not perturb the low-temperature dilute spin-glass state significantly. Because each 
muon on the Y site is surrounded by 6 Mo 5+ ions, the fraction of muons in a non- 
magnetic environment of 0.60(5) implies that at low temperature ( < 10 K) a fraction 
of 0.91(1) of the Mo 5+ sites [ (0.91) 6 = 0.6] are non-magnetic. Hence, only ~ 9(1)% of 
the Mo 5+ carry a net magnetic moment at low temperature. Below 400 mK no further 
increase of the muon relaxation rates is observed down to 50 mK, indicating that thermal 
fluctuations do not play a major role at these temperatures. The spin fluctuation rate 
t and magnitude A can in ideal cases be determined from the field dependence of the 
exponential muon relaxation rate Ai using Redfields formula [30] (see supplementary 
information). The field dependence of the exponential relaxation rate Ai in El at 50 
mK as obtained from the data in Fig. 3b was fitted with Redfields formula (the green 
line in Fig. 3d, see also supplementary material), yielding a fluctuation time scale of 
5(2) ns, with a field distribution A of 2.2(1) mT. This field-strength is as expected for 
muons in El, without directly neighbouring unpaired spins. These muons probe the 
average fluctuating magnetic field from many distant unpaired electrons. 



1 < f3 < 2 (for E2); 



P-M = P mc (t) [/e- Alt + (1 - /) [e~^ P + b 



) 



(i) 
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3. Discussion 

Our analysis of the muon data rules out any major rearrangement of the unpaired spins 
over the available Mo sites within the time scale of the muon experiment of 16 fis. 
Dynamical models of fluctuating valence bonds, where the orphan- or dangling spins 
are mobile, as proposed for SrCr 8 Ga40ig |26j have been explored but can in the case of 
Ba 2 YMo06 be ruled out because of the much larger fraction of dangling spins, while at 
the same time the muon relaxation is more than an order of magnitude slower than for 
SrCrgGa^ig [26]. Such models predict a fast t = muon relaxation for a significant 
fraction of the muons, which is not observed; the t = asymmetry at 50 mK and 120 K 
are comparable. Hence, the low temperature state in Ba2YMo06 is static inhomogenous 
state without any sign of local order parameters down to 50 mK. This state contains 
non-magnetic sites where the neighbouring Mo 5+ Ad 1 electrons must have paired up in 
non-magnetic spin-singlets, and residual unpaired spins. A VBG state, as previously 
proposed [2D], provides a natural explanation for this inhomogeous state. However, a 
larger than previously estimated [2H 120] amount of Mo/Y antisite disorder would also 
explain the presence of residual unpaired spins at fixed locations. Hence, of all data 
so far available on Ba2YMo06 the strongest evidence in favour of a VBG is the heat 
capacity [20] . As pointed out earlier, the pseudogap in the heat capacity [20J below 50 K 
(also visible in the neutron spectra [H]) points to a gradual, strong-glass, freezing at 
temperatures comparable to the magnetic exchange interaction | J\\ of ~ 73 K. Hence, 
it must be the readiness of each 4d electron to randomly form a near-neighbour valence 
bond with one of its 12 near neighbours at high temperature, that prevents freezing 
into a (quantum-) ordered structure at low temperature. This type of behaviour is also 
characteristic of structural glass formers. It is also thought that disorder, even a very 
small amount, can stabilise a VBG state [21] [31], although in other cases disorder is 
known to reduce kinetic constraints towards nucleation of the thermodynamic state. An 
intriguing question that therefore remains is whether the thermodynamic state would 
be a topologically ordered state, in which case the kinetic constraints stabilising the 
VBG are due to topological overprotection, or an orbitally ordered state leading to a 
structural distortion from cubic to tetragonal symmetry at low temperature. It should 
in this respect be noted that the latter is observed in the closely related materials 
Ba2NdMo06 and Ba2EuMo06, while analogues with Gd, Dy, Er, Yb have a magnetic 
susceptibility similar to that of Ba2YMoOg [23]. 

We now turn to the observations on the ~ 8% of residual unpaired spins, that 
account for the low-temperature magnetic susceptibility. The Arrhenius-activated 
behaviour in the low-temperature AC susceptibility indicates that these spins freeze 
into a dilute-spin-glass-like state, with thermally-activated collective dynamics at the 
macroscopic scale. This implies that these spins interact via a weak residual exchange 
coupling (dipolar couplings would lead to freezing at still lower temperatures), and are 
randomly distributed over the Mo sites. It is surprising that even for muons directly 
neighbouring these unpaired spins the relaxation down to 50 mK remains extremely 
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slow; if the residual spins giving rise to the spin-glass response in the AC susceptibility 
were static at the time scale of the muon experiment, this would lead to fast relaxation 
of 2/3rd of the implanted muons well within 1 /xs, as observed in dilute spin-glasses [32J. 
Less extreme examples where these apparently disparate behaviours coexist are known 
amongst geometrically frustrated magnets with ground states sometimes described as 
"quantum disordered" [251 EE]. The slow relaxation of the muons can only in part be 
explained by the large distance between the muon site and the nearest spin moments. 
A full understanding of this common behaviour in low-spin frustrated magnets is still 
outstanding. It is possible that the implanted muons perturb the fragile low-temperature 
spin-glass state, but this raises the question why below ~ 10 K the nuclear spins 
usually remain static on implantation of a muon. Zero-temperature fluctuations might 
also arise due to decoherence that affects ungapped quantum states. As a result the 
quantum glass state is best represented by a density matrix and not as a single wave 
function [12]. This is however unlikely to account for the difference in behaviour at 
microscopic and macroscopic scales. A third explanation that should be considered 
is that the explanation lies exactly in the difference of length scales; that the spin- 
glass freezing observed in the (macroscopic) magnetic susceptibility is truly an emergent 
property, while at microscopic length scales the spin dynamics is dominated by quantum 
fluctuations. These quantum fluctuations are in this case ineffective in equilibrating the 
system to its thermodynamic (quantum-) ordered ground state [I2j HU [13]. A time 
resolved coherent x-ray diffraction experiment on out-of-equilibrium states with local 
order parameters does suggest [33j that quantum fluctuations can at microscopic length 
scales occur at ns time scales, as also suggested by our data. 

4. Conclusion 

Our results confirm that the magnetic low-temperature state in E^YMoOg without 
broken symmetries in the magnetic degrees of freedom is an inhomogeneous state. 
It could be dynamic but at a time scale that is slow compared to the time scale of 
the muon experiment of 10 /xs, with the majority of spins paired into valence bonds. 
That this is a valence bond glass rather than a quantum spin liquid is in this case 
concluded from the electronic heat capacity [20] ■ The /xSR experiment carried out here 
does not distinguish between a spin liquid and a valence bond glass; the relaxation 
is very similar to that reported for many materials presumed to have spin-liquid-like 
ground states [6]. This highlights the importance of the electronic/magnetic heat 
capacity to characterise materials with spin-liquid like ground states, providing for the 
present compound still the strongest confirmation of a VBG, while for the 2D 5 = 1 
antiferromagnet Ba3NiSb20g the magnetic heat capacity clearly points to a spin liquid 
state [9]. In the absence of such data for the well-known kagome antiferromagnet spin 
liquid material herbertsmithite, it remains to be seen whether it is a spin liquid or a 
valence bond glass. The possibility of a VBG due to a small amount of disorder was put 
forward [3T] following neutron spectroscopy data displaying below 120 K temperature 
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independent short-ranged antiferromagnetic dynamic correlations pU], very much like 
the neutron spectrum observed for P^YMoOg [IT] . 
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